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Resume 

Nous nous interessons a la convergence vers sa moyenne spatiale ergodique de la moyenne temporelle d'une 
observable d'un flow hamiltonien a un degre et demi de liberte avec espace des phases mixte. L'analyse est faite 
au travers de revolution de la distribution des moyennes en temps fini d'un ensemble de conditions initiates 
sur la meme composante ergodique. Un exposant caracterisant la vitesse de convergence est defini. Les resultats 
indiquent que pour le systeme considere la convergence evolue en t a , avec a = 0.45 pour alors qu'elle evolue en 
t 1 / 2 lorsque la dynamique est globalement chaotique dans l'espace des phases. De meme une loi a = 1 — /3/2 
reliant cet exposant a a l'exposant caracteristique du deuxieme moment associe aux proprietes de transport P est 
proposee et est verifiee pour les cas consideres. Pour citer cet article : A. Noml, A. Nom2, C. R. Mecanique 333 
(2005). 

Abstract 

Ergodicity, stickiness and anomalous transport We consider the problem of convergence towards spatial 
ergodic average of the time average of an observable defined for a one and half degree of freedom Hamiltonian flow 
with mixed phase space. The analysis is performed by analysing the evolution of the distribution of finite-time 
averages. An exponent characterising the "speed of convergence" is defined. Results indicate that for the considered 
mixed case, the rate of convergence goes as t a , with a = 0.45 while it goes as t 1 ^ 2 when the full phase space is 
chaotic. Moreover a formula linking this characteristic exponent to the one corresponding to transport properties 
f3 is proposed a = 1 — /3/2 and good agreement is found for the considered cases. To cite this article: A. Noml, 
A. Nom2, C. R. Mecanique 333 (2005). 

Key words: Dynamical Systems ; Hamiltonian Chaos ; Anomalous Transport 
Mots-cles : Systemes Dynamiques; Chaos Hamiltonien; Transport Anomal 



Email addresses: Xavier.Leoncini@cpt.univ-mrs.fr (Xavier Leoncini), Cristel.Chandre@cpt.univ-mrs.fr (Cristel 
Chandre), omeziani@yahoo.fr (Ouerdia Ourrad). 



Preprint submitted to Elsevier Science 



May 19, 2008 



Abridged English version In this note we consider the problem of convergence to the ergodic average 
of an observable of a one and a half degree of freedom Hamiltonian flow, and its connection to transport 
properties. The considered observable is the norm of the speed (p, q) defined in phase space, and the 
transport properties of the associated variable, namely the length in phase space of trajectories. 

Assuming ergodicity we conclude that the infinite time limit distribution of time averaged speeds (1), 
is a Dirac. Then when supposing that we have a very chaotic system of Anosov type, we can expect from 
the central limit theorem that the finite time average distribution has a maximum peak which grows to 
infinity as t 1 / 2 . We consider this last property of the height of the maximum of the probability distribution 
function of finite-time averages for a system with mixed phase space and use it to define a characteristic 
exponent a, expecting an algebraic growth as t a . 

In order to link this exponent to transport properties we make the crude approximation of neglecting 
the tails (described in Fig .1) to link the exponent (3 characterising the second moment of transport 
properties with a and obtain the law a — 1 — /3/2, which should be valid for Gaussian transport. 

We then test the ideas numerically with the simple pertubed pendulum (4) for two different values of 
the parameters, one giving us more or less uniform chaos in a given region of phase space, the other one 
presenting a chaotic sea within which regular islands are present (see Fig. 2). Numerical results are made 
using the fifth order optimal symplectic scheme described in [3], using 1024 different initial conditions. 
Trajectories are computed for 10 6 periods with a time step St = T/200. Histograms have been computed 
with a resolution of 5000 bins. The evolution of the averaged speed distribution and the evolution of the 
maximum in log-log plot are represented in Fig. 3. One typically see that for both cases, the exponent a 
can be determined. We find a = 1/2 for the chaotic case while a — 0.43 for the mixed case. More over we 
compute the transport properties for both system and compare the measured f3's with the ones expected 
from Eq. (3). Results are presented in table 1, and a good agreement is found with the extrapolated 
exponents form the a ~ 1 — (3/2 rule. 



1. Introduction 

Dans cette note nous considerons un flot hamiltonien a un degre et demi de liberte H (p, q, t) et perio- 
dique dans le temps et nous nous interessons a la convergence de la moyenne temporelle d'une observable, 
la norme de la vitesse dans l'espace des phases. Plus precisement nous considerons revolution de la 
distribution des moyennes temporelles d'un ensemble de differentes conditions initiales et etudions les 
proprietes d'ergodicite du systeme[l]. Pour une condition initiale donnee i, la vitesse moyenne au bout 
de n periodes s'ecrit : 



On considere alors un ensemble de conditions initiales que Ton suppose equivalentes. Ainsi si la dynamique 
est ergodique chaque Vi(n) va tendre vers la moyenne spatiale quand n — > oo, et done la distribution des 
Vi{n) pour l'ensemble des conditions initiales equivalentes tendra vers une distribution de Dirac. L'objet 
de ce travail est d'etudier comment cette distribution est approchee (si elle existe) et de definir un expo- 
sant caracteristique. Nous considerons dans un premier temps un cas "ideal" d'un systeme Hamiltonien 
globalement chaotique pour lequel les resultats observes sont conformes au theoreme de la limite centrale, 
pour considerer ensuite un systeme ou l'espace des phases est mixte, e'est-a-dire un systeme pour lequel 
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Fig. 1. Estimation de la relation entre l'ecart type et le maximum d'une distribution 

nous avons une mer chaotique peuplee d'ilots de stabilite au sein desquels la dynamique est reguliere. 
Pour ce dernier systeme nous considerons uniquement des conditions initiales dans la mer chaotique. 



2. Definitions 

Nous nous proposons de definir un exposant caracteristique a partir des considerations suivantes. Tout 
d'abord imaginons un systeme chaotique, de type Anosov, tel que Ton puisse appliquer le theoreme de la 
limite centrale. Dans ce cas la variable yi{n) = ^/n,Vi(n) pour n grand va se comporter comme une variable 
aleatoire distribute sur une gaussienne d'ecart type fini. Savoir comment la distribution converge vers la 
distribution limite est l'objet de la theorie des grandes deviations mais cette propriete de convergence 
ergodique implique que la distribution des Vi(n) va tendre vers un Dirac (une gaussienne d'ecart type 
nul). De meme on peut anticiper que l'ecart type de cette distribution va decroitre en l/v^ et que par 
conservation des aires (la probabilite totale est conservee) le maximum de la distribution p ma x(n) va 
croitre en v / "-[2]- Cette derniere propriete permet de definir un exposant 

Pmax{n) ~ n a (2) 

(avec a = 1/2), que nous pouvons a priori generaliser pour des systemes moins chaotiques avec la possi- 
bilite d'avoir a ^ 1/2. 

Interessons nous maintenant aux proprietes de transport du systeme de type Anosov, en mesurant la 
longueur des trajectoires dans l'espace des phases, et done a revolution de la distribution des variables 
Si(n) ~ ^JnTyi(n). Par les memes considerations que celles faites precedemment nous attendons un 
comportement diffusif avec une croissance lineaire de la variance. Maintenant considerons un systeme ou 
le transport est anomal comme par exemple un systeme hamiltonien a un degre et demi de liberte [4] . Le 
transport est dit anomal, dans le sens ou la variance M 2 = ({si(n) — (si(n))) 2 ) ~ wP a un comportement 
non-lineaire. Cette anomalie fait suite aux effets de memoires engendres par le phenomene de collage 
autour d'ilots de stabilite [5,4] qui engendrent une decroissance en loi de puissance des queues de la 
distribution de Sj(n). Dans ce cas en faisant le meme type de raisonnement que celui illustre sur la figure 1, 
un raisonnement negligeant la contribution des queues de distribution et en utilisant la proportionalite 
Vi(n)=Si(n)/nT, nous obtenons 
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Fig. 2. Section de Poincare du hamiltonien (4) pour les parametres suivants : (a) v = 0.5, k = 1, e = 1 et (b) v = 5, k = 2, 
e = 1. 



Ainsi on qualifiera (par abus de langage) le transport de diffusif si a = 1/2, mais aussi on pourra anticiper 
un comportement anomal avec a < 1/2 si le transport est superdiffusif et a > 1/2 si il est sous-diffusif. 
De meme la condition de convergence ergodique a > implique dans le cadre ou la loi (3) s 'applique 
P < 2. 



3. Resultats 

Afin de tester les idees et definitions enoncees auparavant nous considerons le systeme hamiltonien du 
pendule perturbe, qui est egalement le hamiltonien decrivant revolution d'une particule chargee evoluant 
dans le potentiel cree par deux ondes "electrostatiques" dans le referentiel lie a l'une d'elle [6,7,8] : 

2 

H = h A (cos(kiq) + e cos(k 2 q — ojt + ip)) , (4) 
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ou le couple (p, q) designe les variables conjuguees du hamiltonien. Dans la suite et sans perte de generalite 
nous prendrons m — A — k\ = l,ip = et notons k 2 = k et considererons les deux cas (a) v = 0.5, k = 1, 
e = let(b)v = 5,fe = 2,e = l. Nous pouvons noter que sur les differentes sections de Poincare 
representees en Fig. 2 que le cas (a) correspond a une zone chaotique etendue delimitee par des tores de 
rotation reguliers de type KAM, tandis que pour le cas (b) la zone chaotique est peuplee d'ilots au sein 
desquels des trajectoires regulieres sont presentes. 

Nous suivons alors revolution de la distribution des Vi(n), pour cela nous considerons un ensemble de 
conditions initiales prises dans la mer chaotique. L'evolution vers un Dirac de cette distribution pour 
les cas (a) et (b) et revolution du maximum de la distribution sont representees sur la figure 3. A la 
difference du cas (a) nous notons la presence de pics secondaires dans la distribution pour le cas (b). 
Ces pics s'expliquent par le phenomene de collage autour des ilots [4]. Ce phenomene n'est quasiment 
pas observe pour presque toutes les conditions initiales dans le cas (a) (une trajectoire se met a coller 
sur le tore superieur apres 10 5 periodes). Numeriquement, nous avons integre la dynamique en utilisant 
le schema symplectique d'ordre cinq optimal decrit dans [3]. Nous avons considere un ensemble de 1024 
trajectoires differentes que nous avons calculees pendant 10 6 periodes avec un pas de temps St — T/200. 



a = 1-/3/2. 



(3) 




Fig. 3. Haut : Evolution du maximum de la distribution des vitesses moyennes Pmax(n) en fonction du nombre de pe- 
riodes pour les deux cas consideres. L'exposant a caracterisant la "convergence vers l'ergodicite" p m ax(n) ~ n a est mesure 
a a « 0.5~ pour le cas (a) et pour le cas (b) on note deux pentes a f» 0.35 pour n < 2.5 10 3 et a ~ 0.45 pour 
2.5 10 3 < n < 6.5 10 5 , les li gnes en pointille ont une pente 1/2. Bas : A gauche, densite de probabilite des vitesses moyennes 
pour les cas (a) et (b) pour differents temps (t = 160, 320, et 640 periodes). Les distributions ont ete calculees avec une 
resolution de 5000 points a partir de 1024 trajectoires calculees pendant 10 periodes. Le pas de temps est St = T/200. A 
droite, evolution du second moment de la distribution des longueurs parcourues M2, les pentes mesurees sont /3 = 1.55 puis 
(3 = 1.13 aux temps longs. 
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0.5 


0.44 



Tab. 1 

Recapitulatif des exposants observes. 

Nous pouvons noter vu les resultats presentes sur la Fig. 3 que la caracterisation de l'exposant a est 
relativement aisee, les lois etant assez lineaires sur plusieurs ordres de grandeurs. De meme la mesure de 
l'exposant caracteristique du transport associe a Sj, calcule pour les presque toutes 1023 trajectoires pour 
le cas (a) et toutes les trajectoires dans le cas (b) (voir tableau 1) montre que la loi approximative reliant 
les deux exposants est bien verifiee pour les temps long, alors qu'elle ne l'est pas pour des temps courts 
dans le cas ou l'espace des phases est mixte. 

Dans cette note nous avons defini un exposant caracteristique lie a la convergence vers la moyenne 
ergodique d'une observable. Pour les cas consideres, la connaissance de cet exposant permet de predire 
la nature du transport et ce meme pour un cas ou l'espace des phases est mixte. Du point de vue du 
transport anomal nous pensons que cette approche est complementaire de celle liee a la notion de jets 
chaotiques [9] et de e— complexite [10]. 
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